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Abstract 

We explore the idea of accidental dark matter (aDM) stability in the scale invariant local U( 1) b~l 
model, which is a theory for neutrino and at the same time radiatively breaks scale invariance via 
quantum mechanical dynamics in the U(1)b~l sector. A real singlet scalar can be accidental DM 
with an accidental Z 2 , by virtue of both extended symmetries. A U(1)b-l charged complex scalar 
can also be a viable accidental DM due to an accidental (or remanent) Z 3 . They can reproduce 
correct relic density via the annihilations through the conventional Higgs portal or dark Higgs 
portal. The dark Higgs portal scenario is in tension with the LHC bound on Zb~l , and only heavy 
DM of a few TeVs can have correct relic density. In particular, DM may trigger CSI spontaneous 
breaking. The situation is relaxed significantly in the Z 3 case due to the effective semi-annihilation 
mode and then light DM can be accommodated easily. In addition, the Z 3 model can accommodate 
the GeV scale 7 —ray excess from the galactic center via semi-annihilation into pseudo Goldstone 
boson (PGSB). The best fit is achieved at a DM about 52 GeV, with annihilation cross section 
consistent with the thermal relic density. The invisible Higgs branching ratio is negligible because 
the Higgs portal quartic coupling is very small < 10 -3 . 
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I. INTRODUCTION AND MOTIVATION 


Discovery of the standard model (SM) Higgs-like scalar boson at the LHC inspires new 
ideas to understand the following two basic questions about this special member of SM: 
How does it develop vacuum expected value (VEV), namely break the electroweak (EW) 
symmetries? And why does its VEV keep as low as 100 GeV in the presence of quantum 
corrections? Extending the spacetime symmetries of the SM by including classical scale 
invariance (CSI) potentially provides answers to these two questions simultaneously. This 
symmetry implies that the cut-off scale used in the cut-off regularization turns out to be 
a tool instead of a physical scale, and consequently the hierarchy problem may not be a 
physical problem [1], Moreover, the CSI is anomalous and can be broken by undergoing 
dimensional transmutation: via Coleman-Weinberg (CW) mechanism in the perturbative 
region [2-4], or via confining dynamics in the strong coupling region [5]. It opens a chance 
to understand the origin of EW scale from pure quantum mechanical dynamics. In order to 
maintain perturbativitiy of the theory up to high energy scale, it is favored to implement the 
CW mechanism in a hidden sector and then transfer the resulting scale to the SM sector [3]. 
Otherwise, one has to turn to a special structure of dynamics [6]. 

Although not necessary, the hidden sector is supposed to contain a gauge group in order 
to adopt the original CW mechanism. It seems that the local U{1)b-l gauge symmetry 
is a very good case [4], Actually, we do have another strong motivation to extend the SM 
gauge group by introducing U(1)b~l, be., to understand why the observed neutrinos are 
massive but very light. As is well known, incorporating right-handed neutrinos (RHNs) can 
render active neutrinos massive through the elegant canonical seesaw mechanism. In the 
framework of CSI, introducing extra scalar singlets developing VEVs can realize the seesaw 
mechanism maintaining classical scale invariance. Those singlets play an important role 
in CSI spontaneously breaking [7, 8] 1 . Nevertheless, why we need such RHNs and how 
many of them we need may have more profound physical reasons. As is well known, the SM 
fermions are not arbitrary; they are deliberately arranged in order to fulfill gauge anomaly 
cancelations. Similarly, the RHNs, with total number three, gain that kind of legitimacy 
in the local U{1)b~l models. Therefore, the scale invariant B — L model provides a good 
example to address hierarchy problem as well as phenomenological drawbacks of SM in the 
neutrino sector which strongly hint for new physics. 

In this paper we address another motivation for physics beyond the SM (BSM), namely 
lack of a viable dark matter (DM) candidate, in the scale invariant U{1)b~l extension of the 
SM. This topic has been discussed by quite a few groups already, e.g., Refs. [8, 10, 11, 13]. 

1 Scale invariant models that radiatively generate neutrino masses are also studied [9]. 
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In this paper, with the extended spacetime and gauge symmetries at hand, we attempt to 
understand one of the basic questions about DM, namely, why is it stable? The basic idea 
to address this question was proposed in Ref. [ 8 , 11], i.e., DM could be accidentally stable 
due to the extended spacetime symmetries and held content, rather than a protecting ad 
hoc symmetry imposed by hand. This idea of accidental DM (aDM) is motivated by the 
accidental stability of proton in SM. 

In the context of U{1)b~l extension of the SM, it is found that two non-trivial scenarios 
can be realized: 

1. A real singlet scalar is an accidental DM by an accidental Z 2 , by virtue of two extended 
symmetries (classical scale symmetry and the local U{1)b~l gauge symmetry); 

2. A complex scalar charged under U( 1) b-l can also be a viable accidental DM due to an 
accidental (or remanent [14]) Z 3 , which essentially is by virtue of the local U{1)b-l 
but CSI helps to negate other possibilities. 

Both scenarios can produce correct thermal relic DM density via the interactions through the 
Higgs portal to dark Higgs held. But the latter case is in tension with the LHC bound on the 
Zb~l gauge boson from the Drcll-Yan process, and only heavy DM of a few TeVs can have 
correct relic density. Such a heavy scalar DM may lead to an interesting phenomena, i.e., 
spontaneous breaking of CSI could be triggered by DM rather than by Zb-l • The Z 3 case 
is distinguishable from the Z 2 case, since it possesses the characteristic semi-annihilation 
channel, which can relax the tension encountered in the Z 2 dark model with Higgs portal 
scenario. Moreover, a light DM can be accommodated easily in the Z 3 case, as first pointed 
out in Ref. [14], For example, the GeV scale 7 —ray excess from the galaxy center can be 
interpreted by a light DM with mass around 50 GeV and semi-annihilating into a PGSB 
plus DM. 

This paper is organized as follows. In Section II, we study symmetry breaking mechanism 
in the scale invariant U{1)b~l extension of the SM, and its mediation to the SM sector. The 
feature of the Higgs spectra is also discussed. I 11 Section III, we explore accidental DM 
models with Z 2 and Z 3 symmetries and phenomenology therein. Section IV contains the 
discussion and conclusion. Some supplementary materials are casted in the Appendices. 


II. FROM B-L TO ELECTROWEAK SYMMETRY BREAKING 

The minimal classically scale invariant U{1)b~l model (MSIBL) can naturally solve the 
gauge hierarchy problem, by breaking the electroweak (EW) SU(2)l x f/(l)y symmetry 
without introducing the bare mass term for the Higgs doublet. By introducing an extra 
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U{1)b-l symmetry and a scalar field <f> which is a SM singlet but carries U(1)b-l charge 
2 , we can break the local U{1)b-l symmetry through the Coleman-Weinberg (CW) mech¬ 
anism [2]. Besides, the right-handed neutrinos (RHNs) can acquire Majorana masses and 
realize the standard seesaw mechanism as usual after spontanesous breaking of CSI. Then, 
EWSB can be induced by the Higgs portal (~ interaction between the SM Higgs 

field and a new scalar <3> [3, 4], Let us discuss the details of each aspect in the following, 
including the current collider bounds on the new particles. We will focuses on the Zb~l 
gauge boson, which will be highly relevant to the DM phenomenologies. 

A. Symmetry breaking and Higgs spectra 

The MSIBL has two types of new fields, a singlet scalar <3> to break U{1)b-l gauge sym¬ 
metry and three generations of right-handed neutrinos (RHNs) to cancel gauge anomalies. 
The presence of RHNs naturally explains why the SM neutrinos are massive. While U (1) b-l 
breaking naturally furnishes the Majorana mass origins of the RHNs, it moreover fixes the 
charge of the field $ to be +2. With these new particles, the Lagrangian reads 

£ = £ S m - V(H, 4>) - Qa N^NfNi + Y n Jh'N + h.c^j (1) 

We are working in a basis where the mass matrix for the RHNs are diagonal after U(1)b-l 
spontaneously breaking (BLSB). The classical scalar potential with CSI takes the following 
generic form: 

V(H, <L) = ^\H\ 4 - A^|H1 2 |<L| 2 + ^|<L| 4 . (2) 

The scalar potential contains three dimensionless coupling constants. Later we will see that 
among them A h and A h<j> will be almost fixed by two conditions: the weak scale v ~ 246 
GeV and the mass of the SM-like Higgs boson « 125 GeV. Moreover, A will be related 
to gs-L via the hidden CW mechanism in the B — L sector. For the sake of EWSB, the 
coupling constant A^ is assumed to be positive and very small |A^| <C 1. 

Let us detail the above picture of U(1)b-l and electroweak symmetry breaking. It is 
well known that the scale invariant version of the SM fails owing to the heaviness of top 
quark, that leads to a maximum rather than a minimum of the CW potential. To overcome 
this problem, one can consider an alternative scenario where CSI spontaneously breaking 
happens in a hidden sector [3]. This scenario can be realized in a lot of extensions to 
the SM, e.g., the gauge group extension or dark (hidden) sector extensions. In the hidden 
CW approach, some SM singlet in the hidden sector can first develop nonzero VEV via 
the CW mechanism in the hidden sector, and then transfers this scale to the SM sector 
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through the coupling of the singlet scalar to the SM Higgs doublet with a negative coupling 
constant [3, 4]. Obviously, the minimal U(1)b-l model serves as a good example of the 
hidden CW approach 2 . We would like to stress that in case of hierarchical VEVs v -C v^, 
namely when the hidden sector governs the CSI spontaneously breaking, one can turn off 
the cross term A/ l 0 |iL| 2 | < L | 2 and work in the limit of a single scalar field, namely $. This 
approximation is good enough in most cases, except that it may give rise to a discrepancy 
in the Higgs mixing angle. We will come back to this point later. 

Let us first study spontaneous U(1)b-l symmetry breaking by the scalar field <f>. One 
can write <f> = ( 0 c i + + i 4>i) /\/ 2 , where 4> c \ is the classical background, 0 / is the massless 
Goldstone boson (GSB) of BLSB and will be eaten by the U(1)b-l gauge boson (denoted 
as Z' = Zb-l hereafter). While <f>R is the leftover physical particle corresponding to the 
pseudo GSB (PGSB) of CSI spontaneously breaking, it will play a crucial role in the later 
discussion on dark matter. 

There are different ways to demonstrate the dynamics a la the CW mechanism, and here 
we follow the original discussion [2] which yields the particle mass spectra immediately. In 
general, the scalar potential at one loop level can be written as 

r (0 cl ) = Atfa + Bcf a log H, (3) 

where Q is the renormalization scale, which is usually taken to be the minimum (</> c i) = v^ of 
the scalar potential so as to avoid a large logarithm. The logarithmic term resulting in the 
scale anomaly is crucial for triggering CSI spontaneously breaking. In fact, the parameter 
B determines the mass of </>#: 

ml = 8 Bvl. (4) 


The parameter B receives contributions from any fields that acquire masses by coupling to 
the background field </> c i. Specifically to the minimal U{1)b~l model considered in this work, 
the coefficient B is given by 


B 


1 

647r 2 


3 x Q$9b~l 


2 


BW v/2) 4 


(5) 


with Qq, = 2. This expression will be modified if we introduce a new field (e.g., the scalar 
DM considered later) that couples to the $ field. To ensure v$ is a stable minimum point, 
the RHN can not be too heavy: 

B > 0 => j < y/6Q%g 2 B_L o- (6) 


2 Strictly speaking, U{1)b-l is not hidden because all SM matters are charged under it, but we still adopt 
this term in the sense of non-obsecration of Z'. 
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These results are in accord with these of Ref. [4] which adopted the viewpoint of renomal- 
ization group (RG) approach. 

Comments are in order. Firstly, in the above treatment we have turned off the coupling 
A h(j> and reduced the field space into the one-dimensional one. This could be well justified for 
the case of A^ <C 1. Secondly, the parameter A is a combination of A and the loop factors: 

3 x Qlsl-i (-g + InQisLi) - 2 Ajv,i/\/2) 3 4 (-? + In AdJ . 

(V 

In dimensional transmutation, the parameter A is related with the parameter B via the 
minimum condition at v^, namely A = — 2 B, which allows us to express A (at Q = v^) in 
favor of Qb—l and Xn 3 ■ 

A (w) ~ —6Qt a B- L ~ ^ + In QWb-l ^ • ( 8 ) 

Thus, A(u 0 ) ~ ol\_ l <C 1. In case of a very small as-L, which is required for a relatively 
light Z', say around ITeV, A (v^) should be as small as 10” 10 . This may raise doubt on 
the single field approximation in the discussion of CSI spontaneously breaking, because the 
mixing term becomes more important. We will return to it later. Last but not the least, 
being the PGSB of CSI spontaneously breaking, 0-mass is suppressed by loop factors. It is 
illustrative to notice the relation 

_ 3Q%0(B-L f _ 1 12i a2 N,i \ / q x 

m|, 2vr V 6 Q% a%_ L )' U 

where we have defined as~L = Qb-l/^ 77 an< i a N,i = A^j/ 47 r. Ignoring the RHN contribution 
4 , one can make an estimation of the PGSB mass as follows: 

m^, m 0.04(p B _ L /0.1) m Z ' ~ 0.078(p B _ L /0.1) 2 v^. 

It is a few orders of magnitude lighter than Z' and typically lies around the weak scale. But 
in principle this mass can vary within a rather wider region, from a few GeVs to ~TeVs, 
closely related to the ways to fulfill the LHC bounds on Z'. In order to see this, we plot 
several values of m$ in Fig. 1 . 

3 This expression is a little bit different to that in Ref. [4]. The basic reason for this difference is that we 
take a different renormalization scheme. The results will coincide with each other after the redefinition 
of quartic coupling constant X( v <p) A(w^) + Q%a 2 B _ L (l + 6 In Q%g B _ L ). The mass spectrum, which is 
expressed in terms of g B -L, is not affected by this redefinition. 

4 In the following, we shall assume this simplification in order to reduce the number of parameters. 
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FIG. 1: Plots of the LHC and LEP constraints on the parameter space of Z' in the MSIBL. Upper: 
exclusion on the (mz 1 , gB-L) plane; Lower: exclusion on the (v^^b-l) plane. In these plots, the 
two dashed lines show the upper bounds on qb~l , 0.5 and 0.7, respectively. The curves labeled with 
50 GeV, etc., stand for the masses of the PGSB <fi without considering the dark matter contribution 
to CW potential. 

Now we turn our attention to EWSB and consider the change in the Higgs spectrum after 
turning on the cross term. It is a consequence of the second term of Eq. (2) that produces the 
negative mass squared term for H, i.e., m 2 H = ^ v ^ > 0 after CSI spontaneously breaking. 
It then determines the weak scale v to be v = \j2m 2 H /\h as usual. In turn, we have the 
relation 


A h^l = \ h v 2 « m 2 h , (10) 

with rrih the SM-like Higgs boson mass, obtained in the approximation of neglecting the 
mixing effect. This relation helps us to reduce one parameter of the model, A h<j>, for con¬ 
venience. If we adopt the Gildner-Weinberg approach [15] that is suitable for handling the 
CW mechanism in the multi-field space, we would get the same relation via determining the 
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( 11 ) 


so-called flat direction p = pn with n = (n^n^) at tree level: 

hi _n 2 h _\ H _ A 

"} A* A* ’ 


where 0 c i = pn^ and h c \ = with p determined at loop level. Note that the components 
of the PGSB V, just like n, are also determined at tree level. In practice, they are nothing 
but the fractions rth and n^, i.e., V = rihh + n^R. For later convenience, let us define two 
mixing angles 6 and (3 as 


tan 6 


nh 
n ^ 


v 

Vcf, 


1 

tan (3 


( 12 ) 


Then the mixing angle 6 becomes very small, 6 < 0(0.1) for v$ > 3.5 TeV. 

A few comments are in order. Firstly, let us notice that the above expression is different 
from the one derived in Ref. [4] where tan# = cot (3/{ 1 — m^/m^). Note that they coincide 
with each other only in the limit m 2 ^ <C m 2 h with given by Eq. (4). In particular, the 
mixing angle 9 will blow up in their expressions as two masses become degenerate. Since the 
PGSB components should be related to the ratio of two VEVs only, Eq. (12) is a more precise 
expression. Secondly, in this approach, the aforementioned extremely small A problem does 
not bother us because we are now considering a complete potential. Finally, we mention 
that in this model the vacuum stability problem in the SM is almost not affected, because 
the off-diagonal entry of the mass matrix is tiny, and hence the Higgs mixing effect merely 
gives a tiny shift to the diagonal elements. In other words, Eq. (10) works quite well to 
determine A h- In what follows we will denote V as </>. 


B. Z' and (j) at colliders 

The minimal U(1)b-l nrodel predicts two new particles Z' and 0 that can be searched 
for at current/future colliders. Actually, the current LHC bound on Z' is quite stringent, 
because both quarks and leptons are charged under the U(1)b-l gauge group. It has im¬ 
mediate implications for DM phenomenology too. Thus in this subsection we will give a 
up-to-date analysis on the constraints on Z'. We will also briefly comment on the prospect 
on the search for 0. 

Firstly, LEP II measured the cross section for e + e” f f (/ ^ e) above the Z-pole, 

which yielded the following constraint on the heavy Z' [16, 17]: 

m Z 'l9 b-l = Qdg > 6.9 TeV . (13) 

Then one can immediately get a lower bound v$ > 3.5 TeV. Next, the current CMS and 
ATLAS searches for dilepton resonances give a bound on the Z' through the Drell-Yan 



processes (qq —> Z' —* £Jt, with £ = e or /i) [18]. The bound is derived in the ( mz',cr ■ Br) 
plane, with a and Br denoting for a(pp — s- Z') and Br (Z' —s- ££), respectively. We can 
calculate the Z’ production cross section at pp collision and the branching ratios for Z' decays 
into the SM fermions in the MSIBL, and then compare their product to the experimental 
upper bound 5 . In this way, we are able to obtain an exclusion plot in the mz> — Qb-l plane, 
which is of interest in DM phenomenology. From Fig. 1, we have two ways to avoid the 
LHC and LEPII constraints: One way is to assume a lighter Z' and smaller Qb-l] the other 
way is to assume a rather heavy Z' which allows a larger Qb-l■ Obviously, in the former 
case the LHC bound turns out to be more stringent than the LEPII bound, while in the 
latter case the LEPII bound becomes more powerful. Note that imposing the upper bound 
on Qb-l < 0.5 yields a similar bound to the LEPII bound (in the heavy Z' region). We also 
show the LHC and LEPII exclusion plots in the (ry,, Qb-l) plane in Fig. 1. 

The PGSB 0 is another prediction of the MSIBL. But hunting for this new particle at 
colliders is not promising because the mixing angle is very small sin# < 0(0.1). We display 
distributions of 0 in the (m^, sin 6) plane in Fig. 2. can vary in a wide region (in particular 
in the Higgs portal scenario), from GeV to TeV. In the relatively heavier region where sin 6 
is comparatively larger, the most feasible method is to observe a pair of boosted W’s from 
the resonant production of 0: GG — > 0 —> W + W~. For illustration, at the LHC@14 TeV 
and for = 1 TeV and sin# = 0.07, the production cross section of the boosted W + W~ is 
~ O(fb). In the highly decoupling region with very small 6 , 0 can be produced associated 
with Z' via qq —¥ Z'* —)■ Z'<p. However, considering the strigent bound on Z' , the production 
rate should be tiny. 


III. CLASSIFICATION OF ACCIDENTAL DARK MATTER IN MSIBL 

In this section we consider embedding the idea of accidental DM (aDM) in MSIBL model. 
With the aid of CSI and local U(1)b-l symmetry, DM can be realized without much dif¬ 
ficulty. Based on the MSIBL, candidates for aDM can be classified into the following four 
cases: 


• Case A: aDM is a singlet with respect to all gauge groups, and could be a singlet 
scalar S' or a singlet fermion y. However, y must be massless in the MSIBL (which 
offers no singlet scalar with VEV to give mass for y), so we do not consider this case 

5 Note that the Z' involved in LEP II search is off-shell and thus the corresponding bound depends only on 
cross section. This is different from the LHC case, which involves the on-shell Z' and thus is additionally 
sensitive to the branching ratio. Thus the LHC bound can be relaxed by opening new decay channel for 
Z'. 
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FIG. 2: Surviving spectrum. All the points satisfy 0.09 < f lh 2 < 0.12 and are allowed by LUX. 
Top (bottom) left: Z' (0) in the Higgs portal scenario; Top (bottom) right: Z' ((f)) in the dark 
Higgs portal scenario. The perturbativity bounds on gs-L is schematically presented in two cases, 
i.e., the line gB-L = 0.5 and 0.7. 


furthermore. Then S becomes the unique candidate. It conserves an accidental Z 2 
(see Sec. IIIA for detailed discussions). 

• Case B: aDM is a scalar Sx carrying B — L charge X. The trivial case with a peculiar 
A" such that it couples to other fields only via |5x| 2 will be of no interest. Then the 
unique option is X = ±2/3, giving rise to accidental DM (see Sec. IIIB for detailed 
discussions). Fermionic candidates suffer from the gauge anomaly, which probably ask 
for a stack of new particles, a case with much less attraction (the same arguments 
apply to the third and fourth cases below). 

• Case C: aDM is charged under SU(2)i x U(1)y, but electrically neutral. It is the 
neutral component of a multiplet (2 j + 1, Qy ) with j an integer or half integer. As 
shown in Ref. [11], the lower j < 1 cases fail in providing aDM. For example, a triplet 
scalar T = (3,0) allows a term LTL with L the lepton doublet. Although the higher 
j > 2 case is viable, it is the SM gauge symmetries instead of SI and local U (1) b-l that 
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guarantee the accidental Z 2 . Such models actually become the minimal DM scenario 
proposed by Strumia et al. [12]. 

• Case D: aDM comes from a multiplet with double charges, i.e., (2 j + 1 ,Qy,Qb-l)- 
For example, the j = 1/2 case may lead to an inert Higgs doublet. However, unlike 
the Case (B), here we do not have a way to fix the quantum number Qb-l- 

In summary, (A) and (B) give the simplest and the most relevant aDM candidates within 
MSIBL models. In what follows we will study them one by one in two separating subsections. 


IV. REAL SINGLET ADM WITH ACCIDENTAL Z 2 

A real singlet scalar S can be accidentally stable in the MSIBL model, which is explicitly 
seen in its most general interacting Lagrangian 6 7 

Com = ~X. sh S\H^H) - (14) 

where S is odd under a Z 2 symmetry, which is automatic or accidental as a result of the field 
content and symmetries furnished by the model. It is important to note that this feature is 
not broken by quantum corrections. Quantum corrections generate operators which are even 
powers of S and Z 2 is not broken at loop level. It is crucial that RHNs carry U(1)b-l charge 
so that the terms SN? terms violating the accidental Z 2 symmetry are forbidden '. The self 
coupling of S (X s ) is irrelevant to DM phenomenologies but relevant to guarantee that the 
vacuum under consideration is a global minimum, or equivalently dark matter parity Z 2 is 
not spontaneously breaking. For that sake, X s should not be very small (says much smaller 
than 0.1) in the light of the discussion in Appendix B. However, viewing from RGE flow, 
it is in strong tension with a large cross term of the two scalars in Cbm- Eq. (C9) implies 
that if A scj) 1, it will decrease A s fast and make it negative at some scale not far from the 
boundary A. 

The large cross term endangers vacuum stability from another aspect. From the structure 
of 0x so in Eq. (C8), it is clearly seen that, to allow A^u^) as large as 5.0 after running from 
a high scale A TeV, one has to rely on a large gB-L{ A) that can substantially slow the 

6 The model Lagrangian becomes automatically renormalizable due to CSI, if we use the dimensional 
regularization with minimal subtraction. 

7 In Ref. [11] where U(\)b-l is not gauged, and hence an exact accidental Z 2 is problematic after incorpo¬ 
rating RHN to interpret neutrino masses. Other advantages having gauged U(1)b-l are that it does not 
incur large coupling around weak scale, and successful DM phenomenology can be accommodated much 
more readily in the presence of dark Higgs portal S' 2 |<1 , | 2 . 
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reduction of A sc j,(t) from t 0 = 0 to t\ = log(u^,/A). Then A, the quartic coupling of <f>, is 
going to be quickly driven negatively large through the large g%-L ~term in (3\, see Eq. (C 6 ). 
This is not consistent with the dimensional transmutation condition Eq. ( 8 ), which means 
|A(n^)| < 1. In summary, to make the vacuum with only $ acquiring VEV absolutely stable, 
we require either fairly weak couplings or a sufficiently low scale A 8 . In Fig. 3 we show two 
numerical samples of RGE flow with a low (left) and high (right) scale A, respectively. In 
each case, we have tried to get a maximal Xg^v^) by adjusting gs-L and A s . At v<p, it is 
required that A s > 0.1 and |A| <C 1. 



FIG. 3: RGE flows of four couplings with a large A^(u^), left: A = lO 3 ^; right: A = 10 10 u<y The 
perturbativity upper bound at A for the quartic couplings and gB-L are 4-7T and yf&K, respectively. 

Asides from stability, DM models with CSI shed light on understanding another basic 
property of DM, namely the origin of its mass [5, 11, 13]. When specified to the MSIBL, 
aDM mass scale originates from CSI spontaneously breaking in the U{1)b~l sector, just like 
the EW scale. More exactly, aDM actually acquires mass from two dynamical sources: 

m s = \(X sh v 2 + Xsjvl). (15) 

However, in most of the interesting parameter space, practically the second source always 
dominates because of the hierarchy v -C and moreover A s h 1 from the DM direct 

detection bounds. In principle, DM mass can vary in a vast region, from a few GeVs to 
TeVs depending on the quartic coupling A S( j, (and A s h in case of very small X si p). 

Concerning thermal relic density of dark matter, there are two main distinct scenarios 
for this real singlet DM. One is the well known Higgs portal in which A^—term plays the 
crucial role while the A S 0 —term merely provides a (major part) mass of DM. One can give 

8 The most convincible A is Mp\. Then we may say that, in the absence of intermediate thresholds, the 
model is free of hierarchy problem even with gravity. Quantum gravity may leave SI in the IR. 
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the crude condition for the need for Higgs portal. It is obtained simply by requiring the 
annihilation rate of SS —>■ 0</> is short for the typical value, namely a^v < 1 pb. The 
analytical expression will be given later in Eq. (19), from which one can infer that one has 
to fall back on the Higgs portal when < 1.0 x ^ 5 ™ j • By contrast, the other possibility 
is realized in the vanishing X s h limit, for which the single X S( p— term accounts for both DM 
mass generation and DM annihilation (into the PGSB pair) simultaneously. This is a well 
expected scenario since the PGSB can be naturally as light as the weak scale. We dubbed 
this as the dark Higgs portal scenario. But practically this scenario is not dark and it can be 
directly detected sooner or later. Both of these scenarios are simplified, if we decouple the 
two Higgs sectors, but take into account their small mixing, which is compensated by the 
large v makes a difference. It may open new possibilities and we will comment on them. 
In what follows we will go into the details of each scenario respectively. 


A. Higgs portal: an old story 


Although Higgs portal DM has been studied extensively, our arguments to arrive this 
kind of DM model shows more theoretical attraction, as stated before. Here we only brief 
its status facing the bounds from the DM direct detection and Higgs invisible decay, which 
is stringent once kinematically allowed. As for its confronting with the indirect detections 
or global survey, we refer to [19]. In this case, the SM Higgs boson h sm mediates a large 
DM-proton spin-independent (SI) scattering which is given by 


M n ) _ ^sh f I n rn n ( V'' f(n) , 2 f (n)\ 

° SI ~ 4 tt m%mt \ ^ ^ + 9 /tg j 

* n \q=u,d,s / 


(16) 


where [i n = m n ms / (m n + ms) ~ m n is the reduced mass of DM and proton. The values of 
the nucleon parameters etc., can be found in Ref. [20]. In this paper we take the value 
in the bracket, defined as hereafter, to be 0.3 for estimation. 

DM direct detection experiments like XENONIOO and especially the recent LUX [21] 
place an upper bound on the combination A sh/ms- For instance, for a TeV scale heavy 
DM one has crgi < 10“ 8 pb, and hence the upper bound on ms/\ s h > 1-0 TeV. On the 
other hand, this factor enters the annihilation cross section of the TeV scale DM, i.e., av ~ 
(X s h/ms) 2 /32-7T < 10~ 8 GeV~ 2 ~ 0(1) pb, which is just at the correct order of magnitude. 
Thus the heavier DM mass region is still allowed, given a relatively large A s h- But this 
region is expected to be closed in the coming years. In the lighter DM region below m^, by 
virtue of the Higgs resonance, there is a narrow room around m^/2 that is hard to be closed 
by direct detections. Higgs invisible decay yields a loose constraint, if we naively require 
its branching ratio is less than 20% [22], To completely rule out this region one has to rely 
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on indirect detections. We show the results in Fig. 4, which is obtained with the help of 
numerical package Micromega 3.2 [23], scanning the parameter space shown in Table I. 


m z ' 

&B-L 

A sh 

^ S(f> 

[0.5, 5] TeV 

[1 x 10~ 4 ,5 x 10” 3 ] 

[0,1.0] 

[1 x 10 -5 , 6.0] 


TABLE I: Scanned parameters. 



FIG. 4: Status of dark matter under the currently most stringent bound from DM direct search, 
LUX [21], In this plot all the points have good relic density, namely have 0.09 < Qh 2 < 0.12. The 
points in the Higgs portal scenario are around the blue band. While the red points are samples 
from the mixing effect with near 2ms, see details in the text. 


B. Dark Higgs portal: DM as a trigger of CSI spontaneously breaking 

Now we turn our attention to the dark Higgs portal scenario. For simplicity, we turn 
off A s h for qualitative discussion. In the actual numerical analysis we will turn it on and 
scan over it. As stressed before, in most case SS —> (fxf> dominates over other annihilation 
modes. The main contributions to this annihilation come from three Feynman diagrams: 
The contact vertex, the S —channel mediated by 0 and the t/u -channel mediated DM itself. 
The amplitude is given by 

A 2 v 2 

Mss^u — — 2 — 2 SC>> o 2 — 'W- (l”' 7 ) 

— 2 m s 
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We have worked in the non-relativistic limit so that the t— and u —channel contributions 
are equal and independent on the scattering angle. The annihilation cross section times the 
relative velocity is given by 


1 

°U V ~ - x 


1 1 

8n 4 m 2 s 


\M h \ V 1 -r> 




1 X S( p 



(18) 


with = m^/rrig. In the limit <C 1 which works well in this scenario, the final expression 
is simplified to 


a H v « 


1 X S( p 

32^ 


1 rrig 
!6vr vl ' 


(19) 


The second equation was obtained by using an approximation, m| ~ X S( pV^/2. Owing to the 
common origin for DM mass and annihilating dynamics, the cross section is quite sensitive to 
VEV v 0 . A relatively low scale v$ and moreover a larger A sc f, are needed. For instance, even 
if Vtf, is around > 3.5 TeV derived from the LEP II experiments, we still need A s ^ ~ 3.0 in 
order to produce the correct thermal cross section a^v ~ lpb. Recall that low v$ is allowed 
only in the heavier Z' region. Thus the viable dark Higgs portal scenario is supposed to 
dwell in such a slice of parameter space: v$ > 3.5 TeV, gB-L >0.1 and A S( j, > 1.0, which is 
shown in the upper panel of Fig. 6 . It is clear that the LEPII and perturbativity bounds 
could help to rule out a large portion of the parameter space, and in the plot we merely 
show points satisfying v$ > 3.5 TeV and gs-L < 0.7. 

We would like to pause here to discuss perturbativity bounds on two crucial parameters 
gB-L an< 3 A S 0 . At one-loop, the evolution of gauge coupling gs-L is independent on other 
couplings and can be solved analytically, 


gB-i(t) 


1 

(—3«/2tt 2 + 4 9 ^ l( a)) 1/2 ' 


( 20 ) 


Thus, given a sufficiently long RGE course <?b_l(0) < 2.5/\/— 1\. For instance, for A > 1 O 7 U 0 
we have gs-L^v<f>) < 0-6- As for we have shown previously its largeness may render 
vacuum unstable. But if A can be fairly low, says 1000 TeV, A s </>(w) around 5.0 is still 
allowed. In this paper we are not restricted to A = M-p\ and accept a relatively low A, which 
justifies the discussions of this subsection. 

In the heavy DM scenario, an interesting picture may arise, namely, CSI spontaneously 
breaking can be triggered in practice mainly by DM held instead of Z'. This may provide 
a theoretical motivation for introducing a scalar DM held, and one of its immediate conse¬ 
quences is a prediction of a relatively heavy DM. We leave a more general discussion about 
such a scenario in the Section VI. The presence of heavy scalar DM modifies the form of 
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FIG. 5: Dark Higgs portal scenario on the v ^ — qb-l plane. The LEP II excludes the region 
with Vfj, < 3 TeV. The green and red points belong to the scenarios where DM and Z' dominantly 
trigger CSI spontaneously breaking, respectively. We have imposed a perturbativity upper bound 
A< 6.0 in numerical scanning. 


CSI spontaneously breaking. Now, the coefficients B and A given by Eq. (5) and Eq. (7) 
respectively receive new positive contributions from dark matter, 


A B = 


X 2 

A scf> 


647T 2 4 


AA = 


X 2 

A s<j> 


647r 2 4 


- — + In —^ 
2 2 


( 21 ) 


Note that for a complex DM scalar discussed later, there will be an extra factor 2. One 
can estimate when DM contribution to the CW potential becomes dominant over the Z' 
contribution, Eq. (5): 


'W ^ 2y/3Qq,g B -L — 14 9b~l 


( 22 ) 
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FIG. 6: Dark Higgs portal scenario on the — qb-l plane. The LEP II excludes the region with 
va, < 3.5 TeV. The green and red points belong to the scenarios where DM and Z' dominantly 
trigger CSI spontaneously breaking, respectively. We have imposed a perturbativity upper bound 
A s< j) < 6.0 in numerical scanning. 

In this scenario the PGSB mass is strongly correlated with the DM mass, with the ratio 
given by m^/ms ~ -^/A^/47r. Therefore PGSB is typically much lighter than DM, given 
that is constrained by perturbativity. Actually, because A s ^ now is forced to be r\_/ 0(1) 
by correct DM relic density, m^ should be at least a few 100 GeVs. It is different from the 
^-dominant case where m<p heavily depends on gs-L and thus can vary in a fairly wide 
region. This difference is manifested in the bottom panels in Fig. 2. 

In Fig. 6 the green points are those satisfying the condition Eq. (22). They are in the 
relatively large v$ region (> 3.5 TeV), where a heavier my can be realized with a relatively 
small g Z '- Since we have imposed an upper bound A s ^ < 6, v is bounded from above as v$ < 
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6.5 TeV. In our numerical study DM relic density within the region Ddm ^ 2 £ (0.09, 0.12) 
is acceptable. Moreover, we only required that the branching ratio of DM annihilation into 
a pair of 0 is larger than 50% as the definition of dark Higgs portal scenario, and thus 
the actual A S 0 can be mildly beyond the region estimated by assuming the DM-nucleon 
scattering cross section in Eq. (19) to be ~ 1 pb. Then the dark Higgs portal scenario is 
found out to works for A S 0 6 (1.0, 6.0). Interestingly enough, if we take the perturbativity 
bound gs-L < 0-5 (with A = M P[ ) seriously, the red points in the bottom panel of Fig. 6 
are almost excluded. In other words, in the dark Higgs portal scenario DM as a trigger of 
CSI spontaneously breaking is favored by the perturbativity condition. 

Now let us turn our attention to the DM direct detection rate in this scenario, where both 
hsM and (ft can give a sizable contribution to the DM-nucleon scattering. From Eq. (A3) it 
is noticed that the couplings for the S — S — hsM and S — S — (ft vertices are enhanced by A 
and v^, respectively, and the suppression from the Higgs mixing angle can be compensated 
by these factors. More explicitly, we have 


0si 


K4, A 2 m 2 n (An )\ 2 

777,0 J 47 r m|m£ V N J ’ 


(23) 


in the limit A s h A sr f, so that the hsM ~mediation is negligible. The prefactor can be 

traced back to the interference between two channels involving hsM and (ft exchanges in the 
f-channel. I 11 particular, in the degenerate limit ~ m^, the cross section vanishes. By 
contrast, if -C m|, the cross section will instead scale as ~ l/m^, and thus enhanced. 
But in the dark Higgs portal scenario this would never happens, because of the 0-mass 
distribution in the bottom right panel in Fig. (2). Therefore we can neglect this factor and 
then the cross section is enhanced only by A^/A^, compared to that in the Higgs portal 
scenario, Eq. (16). In the large ms limit, the dark Higgs portal scenario is safe in the light 
of current LUX bound. Actually, it amounts to giving the bound: 

1/2 


v<j, 


> 


V 72 A s0 x 1.0 x 


10 - 


a, 


‘pbY 


(n) 
SI,up 


TeV, 


(24) 


which is always satisfied. But there is a good chance to observe or rule out the heavy DM 
in this scenario in the near future. As a matter of fact, from the lower panel of Fig. 6 we 
find that osi is a few 10 -9 pb which is within the sensitivity of the next generation of DM 
direct detection experiments. Besides, that cross section is almost constant. The behavior 
is not surprising because crsi = const x (l — m^/mfy 2 x a^v , where the a^v is the DM 
annihilation cross section that determines thermal DM relic density given in Eq. (19). 

Comments are in order. In the minimal analysis presented in this subsection, we have 
neglected three RNHs. In principle it is possible that (ft couples to RHNs with sizable 
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couplings and hence the annihilating mode SS —> NN dominates over other modes. But 
the qualitative results in our previous analysis are not affected much, since the cross section 
for SS —> NN would scale as ~ \ 2 N m 2 s /v similarly to the previous case Eq. (19) except 
for the enhancement factor X 2 N > 1, which however is offset by a significant phase space 
suppression yi - w/ 


C. Modifications from the mixing effects 


In the previous two sections, concerning the DM pair annihilation, we have considered 
interactions in the individual sectors, respectively the SM and U(1)b-l sectors, ignoring 
the mixing between h and </>. But the h — (f> mixing could make some difference. Actually, 
we have displayed its importance in the DM direct detection. In this subsection, we will 
comment on their effects in DM pair annihilation. 

First of all, from Eq. (A3), let us note that the Higgs portal is modified by the trilinear 
coupling S — S — hsM which receives a contribution from Higgs mixing. The effect can be 
absorbed in the redefinition of A s h —> X s h + Kc/>, but this operation produces an effect in 
the quartic coupling for S — S — hsM — hsM vertex. In the light DM region with ms 100 
GeV, we have A s< j, < C>(10 -3 ), which typically is much smaller than A s h and the resulting 
shift would be negligible. For the heavier DM, the shift may be appreciable and could lead 
to deviations from the standard Higgs portal scenario. For instance, one may arrange a 
cancellation between A s h and A S( j, such that the trilinear term vanishes and only the quartic 
coupling survives. But in our numerical search we do not consider such subtlety. 

Secondly, 0 inherits couplings of hsM and consequently it mediates DM annihilating into 
SM particles. For A s h X^, the two mediators, in the massless limit or in the degenerate 
limit between <f> and hsM, show cancellation in DM annihilation by virtue of orthogonality of 
the rotation matrix from (h, (p) basis to the mass eigenstates. It is in analogy with Eq. (23). 
Obviously, there are two possible ways to avoid this cancellation: 


• The first way is to assume m^ ms so that the 0—mediation is suppressed. Then 
only the hsM —mediation would survive. Nevertheless, such a Higgs portal-like scenario 
is distinct from the conventional Higgs portal scenario, because for A s h <C Awe do 
not have the contact interaction 5' 2 hg M , which is crucial for SS —> hs M hg M to catch 
up with the mode SS — > VV in the limit ms my- In other words, this scenario 
does not respect the equivalence theorem. More explicitly, in the limit of massless 
final states, we have 


a vv v ~ 




1 (A S (j> + A s h) z 


64:71 


mt 


(25) 
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while (Jhh,v is suppressed given A s h -C A S( j> and m| 3> m 2 . Here the parameter 5y 
is given by 8 \y = 2 and 8 z — 1 respectively. Note that the EE mode is enhanced 
by the annihilation into the longitudinal component of E, which produces a factor 
(2m 2 s /m,y) 2 [25]. On the other hand, other modes, including into the SM Higgs boson 
pair, are suppressed by 1/mg and thus are negligible. We rewrite Eq. (25) in a more 
illuminating form: 


cr vv v ~ 


S V 


1 Xgtfi 

32 ^ 


8v (T <fxt> v - 


(26) 


with a 00 V given in Eq. (18). Interestingly, the Goldstone equivalence theorem seems 
to be recovered if we replace /ism —>■ 0. However, this is not true because we have 
dropped the 0—mediation by assuming m$ mg and thus practically a^v = 0. If 
that assumption breaks, then the 0—mediation will spoil the relation Eq. (26). 


• By contrast, the second way is suppressing the hsu —mediation by letting near 2 mg 

so as to resonantly enhance the 0—mediation. As a result, DM mainly pair annihilates 
into a pair of W/Z. The cross section can be calculated in terms of the formula in [24] 


cr vv v ~ 


fn. 1 

mlg mg 



r (h; vv) 


6 v 


fn. 

An m^ 


(27) 


with f r fa 1/(4 — m|/m|) 2 in the narrow width approximation. One can find the 
expression for T(h^ —y VV) similar to that in Ref. [25], taking mass of the virtual 0 
equal to 2 mg. We are working in the limit m| 3> my. As Eq. (26), one can rewrite 
cr V yv ~ 32 8y fgG. Accordingly, when m<f, is not far from the resonant pole, e.g., 
m$ = 1.7mg which gives fn ~ 0.8, the WW channel will become totally dominant and 
can readily have a correct cross section for DM even for X s ^ <0(0.1). As a matter of 
fact, A sr p should be sufficiently small, otherwise mg is so heavy that 0 is not able to 
approach the resonant pole. After all, 0 is a PGSB. The DM-nucleon scattering rate 
is given by Eq. (23). The resonant enhancement in DM annihilating means that the 
factor \ S( j,/m 2 s can be small, and thus ogi may be more or less suppressed. 


We do not make a separated numerical study on the mixing effects, but one can get an 
impression on them from Fig. 4, where the scattered points are due to such effects. For 
illustration, we select the samples in the 0—mediation case, i.e., the case near 2 mg, and 

label them red. As expected, they do have a relatively smaller <tsi. 
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V. U(1)b-l CHARGED COMPLEX ADM WITH ACCIDENTAL Z 3 : SHINING 
LIGHT ON GEY SCALE 


If the DM candidate is a U(1)b-l charged scalar, we can get a trivial aDM with an acci¬ 
dental Z 2 by assigning it a peculiar charge such that it can only couple to <h via |S'x| 2 | < h| 2 , 
as discussed in detail in the previous subsection. However, we can also consider a nontriv¬ 
ial case where the charge of Sx is fixed by coupling to 4>. By virtue of SI, the possible 
interiactions among Sx and $ would be & 3 Sx, & 2 Sx, $S'v|*S'x | 2 and The first term 

evidently renders Sx unstable. The second and third operators (along with |$| 2 $S'x that 
is not listed) actually give identical and Sx up to a conjugate, so again Sx is not stable. 
Therefore, X = —2/3 is the unique option and we get DM with an accidental Z 3 discrete 
symmetry [26]. It survives even after U(1)b-l spontaneously breaking 9 . Now, the relevant 
interaction terms are given by 

$S| + c.c.^ +V(H,$). (28) 

SI and EW symmetries spontaneously break in the same way discussed in Section II A, i.e., 
SI spontaneously breaks mainly at the B — L sector by virtue of the heavy Z B -l and is then 
mediated to the SM sector via the Higgs portal term. In the above Lagrangian, the A 3 cubic 
term of Sx (the subscript will be dropped henceforth) does not contribute to the mass of S. 
So the mass of S again is given by Eq. (15). 

The correct thermal relic density of complex scalar DM S can be realized in this model due 
to newly open annihilation channels involving U(1)b-l gauge interactions and the A 3 —term, 
in addition to the usual Higgs portal integration from term which is present in the real 
singlet scalar DM scenario too. Arguably, among the various contributions to DM pair 
annihilation, that from U(1)b-l gauge interactions are always subdominant to from dark 
Higgs portal, due to the heaviness of Z' or/and smallness of Qb-l • In order to see this, let 
us estimate the order of magnitude of DM annihilating into the SM fermion pairs via Z' in 
the s— channel, 


-£ Zs = X s ^\Sx \ 2 \^\ 2 + Xsh\S x \ 2 \H \ 2 + 


A3 
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a ff v ~ 


E^ 


4 vr ol\ l 
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2„ ,2 




(m|, — 4m |) 2 


Q) 




Y,CfQ 




1 ^s<p J 

32^7 


(29) 


where Vf is the relative velocity between the DM pair in the CM frame when they freeze 
out. The summation is over all U(1)b-l charged fermions / which are lighter than DM. For 
the parameters CfQ 2 , we have C 1 Q 2 = 1, CxQ% = 4 and C q Q 2 q = 1/3 for leptons, RHNs 


9 We also regard it as a remnant subgroup of the local U(1)b-l [14]: U(1)b-l would break down to Z 3 
due to the A 3 term after U(1)b-l breaking from nonzero VEV (<f>). 
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and quarks, respectively. Therefore, this cross section is merely about 10% of a^v which is 
given in Eq. (19). When ms > my, the mode SS* —y Z'Z' is kinematically allowed. But 
its cross section is even smaller than the one given in Eq. (29), and this channel plays no 
important role. 

So, we will concentrate on the A 3 cubic term, which will give rise to a crucial difference 
in DM dynamics, ft opens an effective annihilation channel which is characteristic for Z 3 
models, i.e., the semi-annihilation SS —> S*4> (via the contact interaction) 10 . The cross 
section for this channel is simply given by 


- 


_L^I_ fl 

64n m| V 4 


V5 _ r _± 

2 4 ’ 


(30) 


as usual r$ = ft scales as A|/m|, rather than ~ A§/v| which was the case for the 

previous channel. Therefore, for reasonably large A 3 and light (j) in this dark Higgs portal 
scenario, the lighter DM can easily achieve correct relic density even in the absence of a large 
coupling. So, this scenario survives even for a boundary A = Mp\. It is a distinct feature 
from the dark Higgs portal scenario considered in the accidental Z 2 case, where by contrast 
a quite heavy DM is needed. Obviously, the basic reason is ascribed to the separation of 
DM mass source from the main DM interactions for annihilating. 

We would like to add one comment: annihilation from the A 3 cubic term is not related 
to DM-nucleon scattering at all, even if the Higgs mixing is significant. Of course, here DM 
possesses crsi like Eq. (23) due to the term for mass source, |S'| 2 |$| 2 . But that expression 
should multiply a factor 1/4, because we are considering complex scalar which is non-self 
conjugate. In addition, 0 is favored to be light, so we replace —> m^. Then the estimated 
spin-independent direct detection cross section is given by 


cr S i ~ 0.54 x 10” 11 


m s 

100 GeV 
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/ 5 TeV 

V w 


/ 50 GeV A 

\ m <fi ) 


4 

pb, 


(31) 


which is about two orders of magnitude below the current bound. It is quite sensitive to 
Vfj,. We have taken a larger v<p, which is necessary to get a sufficiently light PGSB with the 
help of smaller gs-L■ This DM may leave its hints in the future detectors, especially for the 
relatively heavy DM near TeV. 

This Z 3 version of dark Higgs portal scenario is distinguished from other scenarios con¬ 
sidered in this paper: it can accommodate the GeV 7 —ray excess in the galaxy center while 
others can not. The excess was claimed five years ago by L. Goodenough and D. Hooper 


10 Actually, there is another process relevant for DM freezing out, namely the three to two annihilation mode 
SSS —► S*S*, which decreases number density of DM also. But it is suppressed by phase space and thus 
is negligible except for a very large coupling of |5j 4 . 
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after analyzing the Fermi-LAT satellite data [27]. Its confidence has been increasing in the 
sequent years [28, 29]. Light DM activities can account for the excess. For example, the 
signal is fit very well by a 31—40 GeV dark matter with a^v = (1.4 —2.0) x 10 -26 cm 3 /s [29]. 
In this paper we have the Higgs portal scenario to accommodate such a DM. However, it 
has been excluded by DM direct detection, see Fig. 4. Alternatively, one can adopt the 
hidden sector approach, where DM annihilates into two on-shell particles that subsequently 
decay into the SM fermions [30, 31]. Here we expect signature from the semi-annihilation 
SS —^ S*4> followed by (p —* bb , similar to the second scenario. In Fig. 7 we display the best 
fit, with x 2 — 30 that corresponds to x 2 /d-o.f. = 1.25, to the photon spectrum. It gives 
ms = 51.0 GeV, = 44.3 GeV and a cross section that is slightly larger than the thermal 
one. The invisible Higgs branching ratio is negligible because the Higgs portal quartic cou¬ 
pling is very small A h<j> < 10 -3 . Actually, Ref. [30] already includes this semi-annihilation 
mode in fitting the spectrum. But there it is merely one of the three annihilating modes 
and thus its role is not as important as in our paper, where the semi-annihilating mode is 
unique. 



1 10 100 


-E 7 /(GeV) 


FIG. 7: The 7 —ray spectrum of best fit point. In this fit we have taken the NFW dark matter 
density profile with 7 = 1.20. Data labeled as diamond stands for the residual spectrum used in our 
fitting. It depends on the diffuse background model, which suffers a large systematic uncertainty. 
Taking it into account leads to different residual spectrum, thus fairly different fit [32], 
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VI. DISCUSSION AND CONCLUSION 


In this paper, we have explored the idea of accidental dark matter stability in the scale 
invariant local U(1)b-l model, which is also a theory for neutrino and at the same time 
radiatively breaks scale invariance via the B — L sector dynamics. A real singlet scalar 
can be accidental DM by an accidental Z 2 , by virtue of both extended symmetries. A 
U(1)b-l charged complex scalar can also be a viable accidental DM due to an accidental 
(or remanant) Z s . They can produce correct relic density via the annihilations from Higgs 
portal or dark Higgs portal, with the latter giving rise to an invisible DM at the direct 
detectors. The dark Higgs portal scenario is in tension with the LHC bound on Zb-l , and 
only heavy DM of a few TeVs can have correct relic density. In particular, it may lead to 
DM triggering CSI spontaneously breaking. The situation is greatly relaxed in the Z 3 case 
which benefits the effective semi-annihilation mode and then light DM can be accommodated 
easily. Additionally, it is able to interpret the GeV gamma ray excess. This paper is based 
on MSIBL, but the main idea can be generalized to any local U(l)x (and nonabelian gauge 
group). Given the weakened constraints on Z', the dark matter phenomenologies can be 
accommodated easily. 

There are some open questions of interest. Firstly, we mention that aDM from the higher 
j > 2 multiplet has an interesting feature, i.e., it must be sufficiently heavy, at least a few 
TeVs, so as to suppress the too fast annihilation via the full EW interactions. Intriguingly, 
this is well consistent with the idea of DM triggering CSI spontaneously breaking, which also 
needs a heavy DM. Secondly, DM with a dark Higgs portal (or more general, dark portal) 
is present in many models. There dark matter, unlike the situation in this paper, may be 
probed only by means of indirect detections. However, the current bounds are far from 
clear. Therefore, it is worthy to employ an up-to-date analysis on this kind of models, with 
the current cosmic ray measurement data at hand. We leave answers to the open questions 
for further investigation. 
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Appendix A: Relevant vertices 


In this appendix we list the relevant vertexes that have been used in the text. Firstly, we 
present the coupling involving 0, except for to dark matter: 


h SM 4> 2 with n H 

/4^sm with 


— 3 log 2 ) B sin 9 — X^ (cot /3 — 2 sin 9) 


— 3 log 2 ) B sin 2 9— 


X H 1 + (2 + 3X h /X H ) cot /3 sin 9 - - sin 2 9 




ytfbff wi t h y<t>f = 


m 


gv<j>VV with gy — —Cy —A sin 9 ( C\y — 2, Cz — 1). 


(Al) 


(A2) 


At the leading order in sin0, the trilinear coupling ^hg M is not affected and still given by 
3A hV. Note that the couplings between two Higgs are derived taking into account the CW 
potential. Numerically, one can replace B with Compared to the A^—term, which 

is ~ sin 3 9, it can be much more significant in particular when ~ O(TeV). In that case 
we have ~ —1.6 Bv < O (10) GeV. In some case it may plays a role. It is safe to make 
the approximation ~ —(1 + 3)A h<pv<p where the numerical factors 1 and 3 originate in 
the cross term |//| 2 |<I>' and quartic term |//| 4 , respectively. Next the trilinear couplings of 
DM to Higgs bosons are given by 


fish = V cos 9(X sh + X S(j> tan 0 tan 9) = v cos 9 ( X sh + A S(j> ), 


fJ’sh C 2 j 

U'SMj 

y^ 2 0, fat = Vcj, cos 9(—X sh cot /3 tan 9 + X S(f> ). 


(A3) 


It is safe to approximate fi S( p ~ Xg^v^cos 9 ~ 2m|/u^, where the second approximation holds 
for DM getting mass mainly from coupling to $, which is true for the heavier DM. As for 
the complex DM case, in the A s h —> 0 limit we hnd the relevant interacting Lagrangian 


£ 


Z3,int 


1 A 


/D 


--j= (cos 00 + sin 9h SM ) + y + c.c. 


+ ^s<!>\Sx\ 2 ^ + l^sh\Sx\ 2 h; 


SM) 


with 


_ a 3 

/D — 


f-lscp X S( pV(p COS 9 , f^sh X S( pV . 


(A4) 

(AS) 


Appendix B: Global minimum in a two scalar system 


It is of importance to check that the vacuum considered in the text is the global minimum 
in the presence of an extra singlet scalar field S. If the vacuum with singlet developing VEV 
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is the global minimum, DM will become unstable. This is a generic problem for the potential 
with multi-scalar. And we demonstrate this in a simplified model with two real singlets Si 
and S 2 respecting Z 2 symmetry, just the case in our paper. Their tree level potential is 

v'(s 1 ,s 2 ) = h S 4 + fi S 4 + hi Sl 2 S2 2. (Bl) 


When one is interested in the case with only S\ radiatively developing VEV triggered by 
a large A12, the potential including radiative corrections can be rewritten as 


V{S 1 ,S 2 ) 


B,S\ ( In ^ — 
' v{ 


+ 


A 


12 C 2 C 2 


^ + Vcw(5 2 )lo=« 


(B2) 


To derive it, we have eliminated Q via the tadpole condition. The effective potential for S 2 
is Vcw(S 2 )\q= V i = A 2 S 2 + B 2 S 2 ln(S%/vl). A 2 and B 12 are as usual functions of couplings 
only, and in the large A12 limit they are given by 


B x « B 2 


A 2 


12 


647T 2 4 


A? 


A 2 


- + - 

24 64vr 2 4 


12 


2 + l n (^i2/2) 


(B3) 


We now show that the desired vacuum with v\ 7^ 0 only is a global minimum with given 
conditions. Simply, its vacuum energy is Ei = —Bivf/2. Consider the other vacuum having 
v 2 7^ 0 only. Its value is related to V\ via the tadpole condition again: 


V2 = v\e x / 4 A2 / 2B2 ' 


(B4) 


Here v\ appears due to the previous choice Q = v\. The corresponding vacuum energy is 
E 2 = —B 2 V 2 / 2 . Therefore, as long as A 2 /2B 2 > —1/4, we will have u 2 < V\ and hence 
E 2 > Ei. For a relatively large quartic coupling X 2 such that A 2 /B 2 > 0, this condition 
is always satisfied. In practice, the minimum of the second term of A 2 is about —0.06 for 
A12 ~ 5.4. Then X 2 > 0.14 makes A 2 definitely positive. Additionally, one can easily see 
that the vacuum with both Si and S , 2 developing VEVs is not a minimum for A12 > 0. In 
conclusion, in this kind of two-scalar system one can choose proper parameters to guarantee 
the vacuum of interest is indeed the global minimum. 


Appendix C: RGEs 

In this appendix we list the relevant RGEs for the gauge and quartic couplings. Moreover, 
we numerically study their evolutions from the low scale up to some high scale A. Requring 
perturbativity if couplings at A has significant impact on the low energy phenomenologies. 
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The RGEs for the gauge couplings are given by the following equations: 


dg Y 

dt 

dgB-L 

dt 

dg 

mix 

dt 


1 41 


16n 2 3 6 
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'9 yi 


167t 2 

1 

167T 2 


32 


41 


12 gB-L "b o 9B-Lgrnix + „ dB-Ldmi 

3 o 

~Q gmix ( 9mix + 2 Qy) + —gB-L ( dmix + 9y) + 12 gB-Ldn 


(Cl) 

(C2) 

(C3) 


with (yfmj-r the kinetic mixing parameter between U(1)b-l and C/(l)y. And it is assumed to 
be zero at the boundary. The RGEs for the quartic couplings are given by the following 
equations: 

d\ h 1 


dt 


167T 2 


12 A ft + 2X hcf) + X sh 
+ X h (12 Y 2 - 9 g 2 - 3 g 2 y - 3g 2 mix )] , 

dXfuj) 1 


l2E t 4 + -g l + -fi'y + -gAr + -g 2 gm, 
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